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Surface exchange and oxygen vacancy diffusion dynamics were studied in double-perovskites 
LnBaCo205,5+5 (LnBCO) single-crystalline thin films (Ln = Er, Pr; —0.5 < 6 < 0.5) by carefully monitoring 
the resistance changes under a switching flow of oxidizing gas (O2) and reducing gas (H2) in the temperature 
range of 250 ~ 800 °C. A giant resistance change AR by three to four orders of magnitude in less than 0.1 s 
was found with a fast oscillation behavior in the resistance change rates in the AR vs. t plots, suggesting that 
the oxygen vacancy exchange diffusion with oxygen/hydrogen atoms in the LnBCO thin films is taking the 
layer by layer oxygen-vacancy-exchange mechanism. The first principles density functional theory 
calculations indicate that hydrogen atoms are present in LnBCO as bound to oxygen forming O-H bonds. 
This unprecedented oscillation phenomenon provides the first direct experimental evidence of the layer by 
layer oxygen vacancy exchange diffusion mechanism. 



Perovskite oxides exhibit a rich variety of interesting and important physical properties such as metal - 
insulator transition, giant magnetoresistance, spin blockade, etc. due to the complex and strongly correlated 
interactions among the charge, spin, orbital, and lattice. Among them, the A- site ordered double perovskite 
cobaltates LnBaCo205.5+§ (LnBCO) (Ln = lanthanide, — 0.5 < 5 < 0.5) have recently attracted substantial 
attention not only due to their unusual electronic and magnetic properties as well as the metal-insulator transition 
but also their high electronic/ionic conductions for a variety of applications in solid oxide fuel cells, gas sensors, 
gas separation and permeation, electrochemical pumping systems, chemical energy storage systems, and many 
others. The dynamics of oxygen surface exchange and the diffusion of oxygen vacancy^ become critical not only in 
governing the novel physical/chemical properties in the strong correlated interaction but also in determining the 
device performance such as energy efficiency, power densities of the batteries and fuel cells^ as well as the 
sensitivity and reliability of sensors. Therefore, it is a critical issue to understand the dynamics of their oxygen 
vacancy behavior both on the surface and in the bulk for multifunctional transition-metal oxides^'^. On the other 
hand, the recent researches on the highly epitaxial thin films of single-crystalline cobalt double-perovskites 
LnBaCo205.5+§ (LnBCO) (Ln = lanthanide, — 0.5 < 5 < 0.5) show various critical physical chemistry properties 
such as ultrafast oxygen diffusivity^ and high sensitivity to chemical environments^"^. Our recent research 
indicates that when the LaBCO films are exposed to a switching flow of oxidizing gas (O2) and reducing gas 
(the mixture of 4% H2 + 96% N2, which will be referred to as H2 for simplicity) in the temperature range of 250 ~ 
800°C, their resistance R changes by three to four orders of magnitude in less than 0.1 s^°. To understand the 
dynamic behavior of oxygen vacancy in the LnBCO systems, we carefully studied the resistance R and resistance 
change AR of thin epitaxial films of single- crystalline LnBCO (Ln = La, Er, Pr) as a function of the flow-time t of 
the oxidizing/reducing gases to find direct experimental evidence that oxygen/hydrogen atoms diffuse ultrafast 
through these films layer by layer via the oxygen-vacancy-exchange mechanism. 

In general, a perovskite oxide ABO3 has a structure in which the BO2 layers alternate with the AO layers. In the 
A-site disordered perovskite LaSrCo206, there occurs only one kind of AO layers, i.e., the Lao.5Sro.5O layers in 
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Figure 1 | Schematic views of the A-site ordered cobalt double peroskite oxides, (a) LaBaCo206. (b) LaBaCo205.5. (c) LaBaCo205. In (a)-(c) the atoms 
are distinguished by colored spheres: Co = cyan sphere, O = red sphere, Ba = blue sphere. La = orange sphere, the half-red/grey spheres in (b) indicate 
the oxygen site occupancy is half, and oxygen vacancy = grey spheres in (c). (d) Sequence of the C0O2, BaO and LnO layers in LnBCO. The 
oxygen-vacancy- exchange diffusion between adjacent LnO and C0O2 layers is expected to be faster than that between adjacent C0O2 and BaO layers 
(indicated by large and small arrows, respectively) because more oxygen vacancies are present in the LnO layers. 



which the La^^ and Sr^^ ions are randomly distributed^\ In the A-site 
ordered double-perovskite LaBaCo206j however, there occur two 
kinds of AO layers (i.e., LaO and BaO) leading to the repeat pattern 
of (Co02)(LaO)(Co02)(BaO) (Fig. la,d). The oxygen- deficient dou- 
ble-perovskites LaBaCo205.5 and LaBaCo205 have structures that 
have oxygen vacancies largely in the LaO layers (Fig. lb,c). The 
oxygen vacancy in Fig. 1 (b) can be formed in the forms of order 
and disorder structures. Recently, highly epitaxial thin films of sin- 
gle-crystalline LnBCO (Ln = La, Er, Pr) were grown on (001) LaAlOa 
by using pulsed laser deposition with a KrF excimer pulsed laser^^"^^. 
The oxidation and redox reactions of these films were monitored by 
measuring their resistance R under a switching flow of O2 and H2 as a 
function of the gas flow time t. The ac conductivity measurements 
indicate that both oxidation under O2 and reduction under H2 can 
occur at the temperature as low as ~200°C. 

Fig. 2 shows the R vs. t plots measured for single crystalline ErBCO 
thin films grown on (001) LaAlOs at various temperatures under the 
switching flow of the reducing (H2) and oxidizing (O2) gases. The 
detailed R vs. t measurements for the reduction and oxidation cycles 
at various temperatures (Fig. 3a,b) show that at a given temperature 
the resistance change of the ErBCO films under oxidation is much 
faster than that under reduction. When the gas flow is switched from 
H2 to O2, the resistance drops down by a few orders of magnitude, 
and this change depends sensitively on the temperature. The resist- 
ance under O2 or H2 is much higher at low temperatures than that at 
high temperatures. Below ~400°C the reduction under H2 occurs 
not in in one step, but in two steps (i.e., a very sharp increase followed 
by a gradual decrease). Under the repetitive switching of the redox 
environments, the reduction and oxidation processes of the ErBCO 
films are highly reversible as can be seen from Fig. 3c,d for the case of 
the redox reactions at 350° C. 

The observed resistance of the LnBCO films under the flow of the 
redox gases can be understood by considering the average Co oxida- 
tion state. The resistance should be low when the average Co oxida- 
tion state is fractional (e.g., +3.5 and +2.5) because it signals the 
presence of mixed valence cobalt ions (e.g., Co^^/Co^^ and Co^^/ 
Co^^, respectively) and hence the occurrence of either hopping with 
low activation energy or metallic conductivity. In contrast, the res- 
istance of the LnBCO films should be high when the average Co 
oxidation state is an integer (e.g., +3 and +2) because electron- 
hopping is difficult in the case of single valence. The average Co 
oxidation state is +3.5 for LnBaCo206, +3 for LnBaCo205.5, and 
+2.5 for LnBaCo205. The Co^^/Co^^ and Co^^/Co^^ can be under- 
stood by doping levels of 0.5 holes and 0.5 electrons per Co^^ ion, 
respectively^^. Our first principles density functional theory calcula- 
tions for LaBaCo206 and various probable structures of LaBaCo205.5 



and LaBaCo205(OH) show that hydrogen atoms are present in 
LnBCO as bound to oxygen forming O-H bonds. Then, the average 
Co oxidation state is +3 for LnBaCo205(OH), and +2.5 for 
LnBaCo204.5(OH), and +2 for LnBaCo203.5(OH)2. Fig. 3d shows 
that, as the gas flow is switched from O2 to H2, the resistance R of 
ErBCO increases to the maximum value and then decreases gradu- 
ally to the equilibrium value. As the gas flow is switched from H2 back 
to O2, the resistance rapidly reaches its maximum value with the 
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Figure 2 | Resistance vs. the gas flow time measured for single crystalline 
ErBCO thin films grown on (001) LaAlOs at various temperatures under 
the switching flow of the reducing (H2) and oxidizing (O2) gases. 
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Figure 3 | Resistance changes vs. the gas flow time, R vs. t, of an epitaxial ErBCO thin film, (a) During the reduction (under H2) cycle at various 
temperatures, (b) During the oxidation (under O2) cycle at various temperatures, (c) At 350°C under the switching flow of the reducing (H2) and 
oxidizing (O2) gases, (d) Zoomed-in view of (c) for one set of reduction oxidation cycles. 



maximum rate of — 10^ Q/s and then quickly decreases to the stable 
value. These resistance changes are attributed to the average Co 
oxidation states on the surface of the films. In pure O2 environment, 
ErBCO can be fully oxidized to become ErBaCo206 with the Co 
average oxidation state of +3.5. As the gas flow changes from O2 
to H2, ErBCO would be reduced to become ErBaCo205(OH) and 
then ErBaCo205.5 (by losing oxygen in terms of H2O) with average 
Co oxidation state of +3, which explains the sharp increase in the 
resistance. A further reduction under H2 would lead ErBCO to 
ErBaCo204.5(OH) and then ErBaCo205 with average Co oxidation 
state of +2.5, which explains the sharp decrease in the resistance. 

At temperatures ranging from 230 to 400°C, the resistance change 
of the LnBCO films as a function of time in each oxidation cycle 
becomes oscillatory. As representative examples, the R vs. t and dR/ 
dt vs. t plots obtained at 260 °C for ErBCO are shown in Fig. 4a,b, and 
those for PrBCO in Fig. 4c,d. The oscillations in the dR/dt curve 
occur in both the Co^^^ Co^^ and the Co^^ Co^^^ oxidation 
steps. Molecular dynamics studies on GdBCO^^ indicate that the 
oxygen diffusion is much easier in the ab -plane than along the c axis. 
The oxygen vacancies in the LnO layers facilitate the oxygen ions to 
diffuse along the c-axis via the oxygen -vacancy- exchange diffusion 
mechanism between adjacent LnO and C0O2 layers and between 
adjacent BaO and C0O2 layers. In general, oxygen vacancies of 
LnBCO are present mainly in the LnO (Ln = La, Gd, Pr, Er, Nd) 
layers rather than in the BaO layers. Thus, in covering every c-axis 
length along the c- direction, the oxygen diffusion will experience two 
different rates as depicted in Fig. Id. The hydrogen diffusion along 
the c-axis of LnBCO would experience a similar situation. Since the 
resistance of an LnBCO (Ln = Er, Pr) film is measured with the 
electrodes attached on the film surface, the measured resistance 
reflects the average Co oxidation state of the surface, which is affected 
by the hydrogen/oxygen diffusion from the surface to the inside of 
the film and by that in the opposite direction. The hydrogen/oxygen 
diffusion along the c-direction oscillates between a faster and a 
slower rate, so the average Co oxidation state on the film surface 
would oscillate hence leading to the oscfllations of the dR/dt vs. t 



plot. The resistance oscillation in LnBCO requires the presence of the 
LnO and BaO layers with different extents of oxygen vacancies, so 
that LaSCO does not exhibit resistance oscillation (see Fig. S2 of SI). 
In the A-site ordered LnBCO, the LnO layer has more oxygen vacan- 
cies than does the BaO layer. Thus the change in the average Co 
oxidation state of a C0O2 layer would be greater when its vacancy- 
exchange involves the LnO layer than the BaO layer hence making 
the diffusion associated with the LnO layer more readily detectable 
than that with the BaO layer by resistance measurements. However, 
when the ionic radius difference between Ln^^ (e.g., Pr^^) and Ba^^ 
becomes smaller, more oxygen vacancies may exist in the BaO layer, 
so that the diffusion associated with the BaO layer can become vis- 
ible. It is noted that the dR/dt oscfllation of ErCBO is slightly differ- 
ent from that of PrBCO. Each dR/dt osciflation of PrBCO has two 
components (Fig. 4d), but that of ErBCO has mainly one (Fig. 4b). 
The larger oscillation peak of PrBCO is related to the diffusion 
through the PrO layer, and the smaller one to that through the 
BaO layer. The oscillation peak of ErBCO is related to the diffusion 
through the ErO layer, and the diffusion through its BaO layer is not 
observed most probably because the extent of its oxygen vacancy is 
very small. The atomic- exchange diffusion has been found for several 
pure metal surfaces such as Pt, Ir, and AP^"^^. However, there has 
been no direct experimental evidence for the vacancy- exchange dif- 
fusion taking place through the bulk of a crystalline material 
although this mechanism is often invoked to describe oxygen dif- 
fusion in various materials^^. The resistance oscillation observed for 
LnBCO (Ln = Er, Pr) is the first experimental evidence for the 
occurrence of the oxygen-vacancy-exchange diffusion mechanism 
through the bulk of crystalline oxides. 

The diffusion of hydrogen or oxygen along the c-direction through 
the lattice of ErBCO can be analyzed by Pick's second law^^ for the 
one dimensional diffusion by considering a thin layer of thickness L 
eqmlibrated under a partial pressure Pi of redox gas. A sudden 
change of the partial pressure to another value P2 wfll cause the redox 
gas to diffuse eventually leading to the final state. This relaxation 
process can be monitored by resistance measurements. Given the 



SCIENTIFIC REPORTS | 4 : 4726 | DOI: 1 0.1 038/srep04726 



3 



C] 1.8x10' 





11600 



11605 



11610 



11615 



11620 



73985 73990 73995 74000 74005 74010 



Time (s) 



Figure 4 | (a, b) R vs. t and dR/dt vs. t plots for the epitaxial ErBCO films, and (c, d) those for the epitaxial PrBCO thin films, taken for an oxidation cycle 
at260°C. The insets in (a) and (c) show the zoomed-in view of the R vs. t plot in the regions oft = 11603- 11608 and 7392 -7412, respectively. The inset of 
(b) shows the Arrhenius plot, log D vs. 1000/T, based on the data of Table 1, with the negative and positive slopes indicated by solid and 
dashed lines, respectively. 



carrier densities at the initial and final states as Ci and C2, respectively, 
the carrier density c(x, t) at the distance x from the center of the layer 
at the time t is related to the diffusion coefficient D of the redox gas 
(see SI). If the conductivity is dominated by one type of charge 
carrier, the mean conductance a^Cx, t) is approximated as 



02 — Oi 



-3/2 



where ai and (J2 are the conductance at the initial and final states, 
respectively, and this expression is valid when t is small compared to 

the relaxation time t = — — . Then, the conductance change Aam 
during a short time interval At is given by 



(2) 



where 0 = 211 ^^^{02 — g\)\/\/t. We analyze the resistance oscilla- 
tions observed during the oxidation cycle at various temperatures 




ranging from 240 to 400°C by employing Eq. (2) for each oscillation 
peak to deduce the associated relaxation times x, and subsequently 
the diffusion coefficients D by taking L = 2c where c is the c-axis 
lattice constant (see SI). Table 1 summarizes our results obtained at 
various temperatures ranging from 240 to 400 °C. 

The dR/dt oscillations in ErBCO are absent above 320 °C for the 
Qq2.5+ _^ Q^3+ oxidation process, and also below 280 °C for the Co^^ 
Co^^^ oxidation process (Table 1). This indicates that the Co^^^ 
Co^^ oxidation involves largely hydrogen diffusion, and this dif- 
fusion becomes too fast above 320 °C to cause dR/dt oscillations, and 
that the Co^^ Co^^^ oxidation involves largely oxygen diffusion, 
which becomes too slow below 280°C to cause dR/dt oscillations. 
This interpretation is consistent with our observations that the dif- 
fusion coefficients D for the Co^^^ Co^^ oxidation (—10"^ cmVs) 
are much greater than those for the Co^^ Co^^^ oxidation 
(—10"^- 10"^ cmVs), and the temperature -dependence of D is quite 
weak in the former oxidation process, but it is quite strong in the 
latter oxidation process. Thus our data show that the hydrogen dif- 
fusion is several orders of magnitude faster than the oxygen diffusion 



Table 1 | The relaxation times i and the diffusion coefficients D deduced from the dR/ dt curves in the oxidation cycle at various temperatures 

T(°C) t(xIO-^^s) D(x10-2cmVs) t(x10-i°s) D(xlO-^cmVs) 



240 6.56 ±1.86 1 .64 ± 0.46 

260 2.44 ±0.51 0.62 ±0.17 

280 5.49 ± 0.84 1 .37 ± 0.24 1 .04 ± 0.27 1.45 ±0.19 

300 1.09 ± 0.07 1.39 ±0.08 18.6 ± 2.3 0.082 ±0.013 

320 4.92 ± 0.69 0.31 ± 0.05 49.5 ± 0.2 0.031 ± 0.000 

350 - - 3.52 ± 0.02 0.430 ± 0.003 

400 - - (9.25 ± 0.07) x 1 0"^ 1 63 ± 2 
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in ErBCO. The oxygen diffusion in ErBCO is faster than that found 
for other oxides by several orders of magnitude under similar tem- 
peratures (see SI). The oxygen diffusion in ErBCO is faster than that 
in YSZ, Gd:Ce02 and other similar perovskite oxides by several 
orders of magnitude^^"^^ while the hydrogen diffusion in ErBCO is 
comparable in rate to the silver diffusion in a-Ag2+5S and a-Ag2Te^^. 

In thermally activated electrical conduction, a is related to the 
activation energy Ea as a oc exp(-Ea/RT). Since a is proportional 
to D, we plot log D vs. 1000/T using the data of Table 1 to determine 
Ea (see the inset of Fig. 4b). During the Co^^^ Co^^ oxidation 
process, the slope of the log D vs. 1000/T plot is slightly positive. This 
indicates that the ErBCO films behave as weakly metallic, which is 
consistent with our interpretation that the Co^^^ Co^^ oxidation 
largely involves hydrogen diffusion. Thus the activation energy for 
hydrogen hopping would be small (compared with that for oxygen 
hopping). During the Co^^ Co^^^ oxidation process, the log D vs. 
1000/T plot changes from a negative slope above 280°C yielding Ea = 
1.25 eV to a substantially positive slope below 280°C. This activation 
energy is much lower than that found for other oxide materials with 
disordered oxygen vacancies. 

In summary, under the switching flow of O2 and H2 gases, the 
resistance change measured for epitaxial thin films of single-crystal- 
line LnBCO (Ln = Er, Pr) as a function of the gas flow time t exhibit 
oscillations during the oxidation cycle under O2. This manifests that 
the ultrafast layer-by-layer exchange diffusion of O2 and H2 in 
LnBCO takes place by the oxygen- vacancy- exchange diffusion. 
Especially, the hydrogen ions can superfast diffuse in the ordered 
LnBCO systems may pave a new way for the studies of hydrogen in 
the ordered vacancy systems. These results suggest that the LnBCO 
systems are an excellent candidate for the development of low or 
intermediate temperature energy conversion devices. 
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